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ABSTRACT 

We investigate the environment dependence of dark matter halos in the symmetron modified gravity 
scenario. The symmetron is one of three known mechanisms for screening a fifth-force and thereby 
recovering General Relativity (GR) in dense environments. The effectiveness of the screening depends 
on both the mass of the object and the environment it lies in. Using high-resolution N-body simulations 
we find a significant difference, which depends on the halo's mass and environment, between the lensing 
and dynamical masses of dark matter halos similar to the /(i?) modified gravity. The symmetron can 
however yield stronger signatures due to a freedom in the strength of coupling to matter. 



1. INTRODUCTION 

Many theories of high energy physics, like string the- 
ory and supergravity, pr edict light gravitationally cou- 
pled scalar fields (see e.g. Binetruy ( 2006 ) ; Linde (2008j ) 
and references therein), ihese scalars may play the role 
of dark energy (quintessence) . If these scalar fields have 
non-minimal coupling to matter fields, then they could 
mediate extra forces and are potentially detectable in lo- 
cal experiments and from observations on cosmological 
scales. 

Several laboratory and solar system experiments over 
the last decades have tried to detect a sign of such fun- 
damental coupled scalar fields I Adelberger|2002[ Hoskins 
120071 Bertotti et al.||20U3p, but 



et al. 


1985 


Decca et al. 


2007 



the results so far show no signature of therm Naively, 
the results of these experiments have ruled out any such 
scalar fields, which also can have an effect on the large- 
scale structure of the Universe, unless there is some 
mechanism which suppresses the scalar fifth-force on 
small scales where the experiments are performed. One 
should keep in mind that GR is only well tested on length 
scales ranging from millimeters to the size of the solar 
system. Comparing this to the size of the horizon, this 
leaves a wide range of scales where there could be devi- 
ations from GR. 
To this day we know three su ch types of theoretical 



mechanisms (see Khoury (20101 for a review) that can 
explain why light scalars, it they exist, may not be visi- 
ble in experiments performed near the Earth. One such 
class, the chameleon mechanism dKhoury fc Weltman 



20041 IBrax et al.|[2004l [C lifton ct al. 200 5| |Mota k Mar- 
row||2004a|b[ [Brax et al.]|2010bf |Mota fc Shaw|2007p , op- 



erates when the scalars are coupled to matter in such a 
way that their effective mass depends on the local matter 
density. In regions where the local mass density is low, 
the scalars would be light and deviations from GR would 
be observed. But near the Earth, where experiments are 
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performed, the local mass density is high and the scalar 
field would acquire a heavy mass making the interactions 
short range and therefore unobservable. This mechanism 
is the reason why f{R) modified gravity can lead to vi- 
able cosmologies and still evade local gravity constraints. 

The seco nd mechanism, the Vainshtein mechanism 
(|yainshtem| |1972[ [Deffayet et al.||20Q2[ |Arkani-I Iamed 



jef al._2003j ), operates when the scalar has derivative self 
couplings which become important near matter sources 
such as the Earth. The strong coupling near sources 
essentially cranks up the kinetic terms, which translates 
into a weakened matter coupling. Thus the scalar screens 
itself and becomes invisible to experiments. This mech- 
anism is central to the phenomenological viability of 



theories (Dvali et al.||2000| de Rham et al. 


2008 


Nicolis 


et al.||2008; Hinterbichler et al.ii2010| |Mota et al 




Gabadadze 


2009; dc Rham 2010; Brax et al.||2011ap. 



is the symmetron mechanism (Hinterbichler fc Khour; 



20T0l [H interbichler erar]|201ll |01ive fc Pospelov] 



Brax et al.(|2011bj jCannouji et al.[[2010p. 
nism, the vacuum expectation value (VE 



loury 



In this mecha- 
V) of the scalar 
depends on the local matter density, becoming large 
in regions of low mass density, and small in regions of 
high mass density. The scalar couples with gravitational 
strength in regions of low density, but is decoupled and 
screened in regions of high density. This is achieved 
through the interplay of a symmetry breaking potential 
and a universal quadratic coupling to matter. A simi- 
lar screening mechanism applies for the envi ronmentally 



dependent dilaton model Brax et al. (2010a 



In vacuum, the scalar acquires a nonzero VEV which 
spontaneously breaks the Z2 symmetry <j) ^ —<j). In the 
regions of sufficiently high matter density, the field is 
confined near = 0, and the symmetry is restored. The 
fifth force arising from the matter coupling is propor- 
tional to (/), making the effects of the scalar small in high 
density regions. Because of this effect, dark matter halos 
in high-density regions will produce different scalar fifth 
forces compared to those in low-density regions. 
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This effect has been studi ed for the case of .f(R) gravity 
(Cham eleon mechanism) in Zhao et al. (2011b); Schmidt 
|2010[ ) and the Dvali-Gab adadze-Porratti (UGP) model 
( Vainshtein mechanism) in Schmidt] (2010). It was found 
that in the DGP model the screening ot iialos is almost 
independent of environment while in /(i?) gravity there 
can be a significant environmental dependence. Another 
signature that has been found recently in f{R) an d sym- 
metrqn mod els is that the luminosity of galaxies (Davis 
et al.|[2011b ) might also depend on the environment 



Kecent work on the symmetr on has focused on back- 
ground cosmology, linear Brax et al. ( 2011b|) and non- 
linear structure formation j Pavi^et al. (pUllal, and also 
made halo scale predictions Clampitt et al. (2011). In 
th is paper we use the h igh-resolution IN -body simulations 
of Davis et al. (2011a) to study the environment depen- 
dence ot dark matter halos in the symmetron modifie d 
javity scenario (see for instance IZhao et al.l (I2011a| 
Li & Zhao (2009, 2010); Li et al. (2011); Ferraro ct a 
011^;|(Jyaizu ct a l.^^2008J; Schmidt j ^2009| ) ; |Hellwing fc 



uszkiewicz ( 2009 ) ;lBrax et al. (2011c|Tor N-body 



sim- 

ulations within other scenarios of modified gravity). 

This paper is organised as follows. In Section ([2]) we 
recall the main properties of the symmetron model which 
are relevant for our analysis. Then, in Section ([3]) we in- 
troduce the dynamical and lensing masses of a nalo and 
explain how these are obtained from the N-body sim- 
ulations. In Section (|4| we explain how we define the 
environment of a halo m our analysis. The main results 
are shown and discussed in Section ([5| , and we also com- 
pare our simulation results to semi-anaytical predictions 
in Section ([6]). Finally we summarize and give our con- 
clusions in Section (Q. 

2. SYMMETRON REVIEW 

The symmetron modified gravity is a scalar field theory 
specified by the following action 




pi-^W)'-m 



5* = 



(1) 

where g is the determinant of the metric (7^^, R is 
the Ricci scalar, ipi are the different matter fields and 
Mpi = I/VSttG where G is the bare gravitational con- 
stant. The matter fields couple universally to the Jordan 
frame metric g^jy, which is a conformal rescaling of the 
Einstein frame metric g^^i, given by 



(2) 



The equation of motion for the symmetron field (f> follow- 
ing from the action Eq. ([T]) reads 



(3) 



where the potential is chosen to be of the symmetry 
breaking form 



(4) 



and the coupling is quadratic in • 
the (/) —(j) symmetry 



to be compatible with 



A(0) 



M 



(5) 



The effective potential can then be written as 



(6) 



from which the range of the scalar field (i.e., the range 
of the resulting fifth-force) can be found as 



A. 



(7) 



The range of the field in vacuum, denoted Aq, is given by 
Ao = 

In high-density regions where p„ > p^'NP the effective 
potential has a minimum at = 0. The fifth-force, given 
by 



= VA(0) = ^V0, 



(8) 



is proportional to 4> and will be suppressed. In vacuum, 
or in large underdensities, where pm ^ fi^M^ the (j) — >■ 
—(f) symmetry is broken and the field settles at one of the 
two minima (f> — (f>Q = ±/i/\/A. The fifth- force between 
two small test masses in such a region will achieve its 
maximum value compared to gravity. 



Fn 



2M, 



pi 



dlnA 



(9) 



4>=(t>o 



For very large bodies in the sense that 

2 



a 1 ^ 2^ ( 1 « 1, (10) 



Pbody \ -^body 



the s ituation is quite different ( Hinterbichler & Khoury 
2010 1. Here the symmetry is restored in the interior of 
the body and the fifth force on a test mass outside be- 
comes suppressed by a factor a~^: 



(11) 



If the body lies in a high density environment where (f) — 
0Gnv < (l>o the fifth-force will be further suppressed by a 
factor {(f>em / 4'of' ■ Thus there are two ways a body can be 
screened: it can be large enough as to make a"^ <C 1 or 
it can be located in a high density region where (f)cnv ^ 
00- The latter in particular leads to an environment 
dependence of the fifth-force in a dark matter halo which 
we will investigate in the next section. 

Instead of working with /i, M and A we chose to define 
three more physically intuitive parameters L, f3 and zgsB 
which are the (vacuum) range of the field in Mpc/h, the 
coupling strength to matter and the cosmological redshift 
where symmetry breaking takes place on the background 
level respectively. The conversion to the original model 
parameters are given by 



L ^ 



Aq 
Mpc/h 
(f)oMpi 
M2 

- ^^ssb)^ 



3000i7o 

PmO 



(12) 
(13) 
(14) 



3 



3. THE DYNAMICAL AND LENSING MASSES 

In any universally coupled scalar-field theory, like the 
symmetron, we have the choice of describing the dynam- 
ics of the model in two mathematically equivalent frames 
defined by choosing either or g^^, in Eq. ^ as the 
space-time metric. In the Einstein-frame, the one de- 
scribed by Eq. ([!]), gravity is described by standard GR, 
but the geodesic equation is modified compared to GR: 



- ap 



X 



d4> 



^'^^ + 2(j)^pxf^x^') . (15) 



In the Jordan-frame gravity is described by a scalar- 
tensor-like modified gravity theory, but the matter par- 
ticles follow the geodesies of the space-time metric g^i, 



0, 



(16) 



where T is the Levi-Civita connection of g^^. The pre- 
dictions of the theory is usually easier to derive in the 
Einstein-frame and the corresponding quantities can be 
found in the Jordan-frame by performing the transfor- 
mation Eq. ([2]). 

Working in the Einstein-frame and the conformal New- 
tonian gauge, the line-element can be written as 

ds^ = a^{ri) [~dri^{l + 2<^n) + (1 - 2<^N)dx^] , (17) 

where $Ar is the usual Newtonian potential. Transform- 
ing to the Jordan-frame using Eq. ([2]) we find 

ds^ = 0^(77) [-dry2(l + 2$) + (1 - 2^)dx^] , (18) 

where 



$ ~ $Ar + M((?!)), 



with 



^ A{c^) - 1 



M 



(19) 
(20) 

(21) 



Note that we have neglected a ternj^ 2$Ar(5A(0) ^ <i>Ar 
in the equations above. In the solar-system deviations 
from GR are often phrased in terms of the so-called PPN 
parameter 7. In the case of the symmetron we have 



N 



N 



5A(j)) 



= 1 - 



25A{4>) 



(22) 



The sol ar system constraints for the sym metron were de- 
rived in Hinterbichler & Khoury (20101 and gives a con- 
straint on the range of the field a nd the symmetry br eak- 
ing redshift: L(l + zssb)^ < 2.3 pavis et al.|2011a[ ). 

Typically observations of e.g. clusters probes forces 
(gradients of the potentials) instead of the potentials 
themselves and different observables are related to dif- 
ferent combinations of the potentials. The fifth-force po- 
tential is given by the difference in the above two poten- 
tials 



$ - * 
= ^— = M(0). 



(23) 



^ For the symmetron it was shown in [Davis et al] ( |2011a[ | that 
SA{4>) < SA((l>o) ~ which for the values of /3 < 0(1) we 

are interested in is much less that one. 



Lensing on the other hand is affected by the lensing po- 
tential 

(f> + 

$+ = = (24) 

which satisfies the Poisson equation 

V2$+ ^ AnGa^Spm- (25) 

This is the same equation as in GR since the action of the 
electromagnetic field is conformally invariant and thus 
photons do not feel the scalar fifth-force. In general, 
there will also be a contribution from the clustering of 
the scalar-field AirGa^SV {(/)), but in our case this term is 
negligible as the difference in the clustered and unclus- 
tered energy density of the scalar field is always much 
less than the energy density of matter in a halcj^ We 
define the lensing mass as 



Mr - 



1 

47rGa2 



(26) 



which is the actual mass of the halo. It is determined 
form the N-body simulations by counting the number of 
particles within a given radius. For spherical symmetry 
we can use Stokes theorem, JV'^^+dV = / V$+ • dS = 
47rr^ ^^^^ , which gives 



M^lr) oc r' 



dr 



(27) 



The dynamical mass MD{r) of a halo is defined as the 
mass contained within a radius r as inferred from the 
gravitational potential i.e. 



MD{r) = 



1 



47rGa2 



(28) 



where the integration is over the volume of the body out 
to radius r. For spherical symmetry we can again use 
Stokes theorem on the r.h.s to find 



M]j{r) (X r 



2d$(r) 



N 



dr 



\ dr AP dr 



(29) 



The terms in the brackets are recognized as the sum 
of the gravitational force and the fifth-force. In our 
N-body simulations we measure Mo of a halo by first 
using a halo-finder to locate the particles which make 
up the halo and bin them according to radius. Then, 
we calculate the average total force in each radial bin 
by summing over all the particles in the bin. Note that 
the force obtained in this way can have a contribution 
from the particles outside the halo. For spherical 
symmetry, this contribution largely cancels out and 
we are left with the total force produced by the halo 
itself. Observationally, Mo can be determined from 
measure ments of e.g. th e velocity dispersion of galaxies 
in halos ( |Schmidt|[20T0l ). 



In GR the lensing mass is the same as the dynam- 
ical mass, but they can be significantly different in 

The potential energy of the scalar field satisfies |<5y(0)| < 



1 2SSB S; this term is negligible compared 



|y(0o) - vm -- 

2SSb)^10-6. For , 

to the energy density of matter in a halo. 
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modified gravity. We follow Zhao et al. (2011b) and 
define the relative difference 



AM(r) 



Ml 



1 



d^-/dr 
d^+/dr' 



(30) 



This allows us to quantity the difference between the two 
masses in the simulations. In GR we have Am = while 
in the symmetron model Am will vary depending on the 
mass of the halo and the environment it lies in. The 
theoretical maximum is achieved for small objects in a 
low-density environment where the screening is negligible 
and reads (see Eq. l9| 



A5;^-(r) = 2/32 



(31) 

In the next section we will apply these to the N-body 
simulation results. 

4. DEFINING THE ENVIRONMENT 

The environment a halo lies in can have large effects on 
the fifth-force that operates by the halo since the fifth- 
force is directly proportional to the spati al g radient of the 
square of the local field value (see Eq. (29)). This value 
is small in high-density regions and this will provide the 
halo with an additional screening to the self-screening 
due to its size or mass. 

As a result, when looking for an environmental depen- 
dence it is crucial to choose a definition of the environ- 
ment that does not correlate heavily with the halo mass. 
The quantity one chooses should also allow for an easy 
determination both in our simulations and in observa- 



tions. Such a quantity was found in [Haas et al. (|2011| 
and used in the same an alysis as we n 



rave d one, but tor 

the case of f{R) gravity ( Zhao et al.||2011b ). This quan- 
tity, 

Unj = , (32) 

is defined as the distance to the A^'th nearest neighbor 
whose mass exceeds / times the halo under consideration 
divided by the virial radius of the neighboring halo. A 
large value of D indicates that the halo lives in a low- 
density environment in the sen se that it has no larger 
halos close by. It was shown in Haas et al. (2011) that 
the quantity D = Di i represents the local density well 
and is almost uncorrelated with the halo mass. We have 
explicitly checked that this is also the case for our simu- 
lations. 



We follow Zhao et al. (2011b) and define a high den- 
sity environment as logj^Q U < 1 and a low density en- 
vironment as logj^Q -D > 1. Halos in a low- and high- 
density environments will be called isolated halos and 
clustered halos respectively. In order to study the varia- 
tion of Am with halo mass wc will say that a halo with 
Mi/(Msun/^) > 5 X 10^2 is a large halo, while a halo 
with 5 X 10" < ML/{Msnn/h) < 2 X 10^^ is a small halo. 
The lower limit comes from the fact that smaller halos 
in our simulations are not well resolved (less than ~ 500 
particles per halo) and will therefore not be used in this 
analysis. This choice, arbitrary as it might seem, is made 
so that we have approximately equal numbers of halos in 
each of the two categories. 

In Fig. m we show the 3D halo lensing mass distribu- 
tion together with the corresponding value of A m in the 



simulation box for one of the simulations. Each sphere 
represents a halo; in the left panel the size of the spheres 
is proportional to the halo lensing mass Ml and in the 
other two panels it is proportional to Am; the color indi- 
cates the environment of the halos as illustrated by the 
legend. This plot shows the environmental dependence 
clearly: in high density environments the value of Am is 
generally smaller than in low density environments, and 
the definition of the environment used here is capable of 
capturing this behavior fairly well. 

5. RESULTS 



We use the N-body simulations of Davis et al. ( 2011a[) 
which have been performed using a modified version (Ei 



2011) of the pub hcly available N-body code 
nebe et al.||2001 ). The simulation suite con 



&; Barrow 
MLAPM 

sists of six simulations with different model parameters 
{/?, L, zssb} shown in Table (fTl), and we have calculated 
Am for the halos found in tncse simulations. The ini- 
tial conditions for the symmetron models are the same 
for each simulation and allows for a direct comparison 
of the effects of the different parameters in the theory. 
The halos in the simulation have been found with the 



Mvlr 



halo finder MHF ( |Gin et al.||2004 [) using the definition 
M(r34o) 



M340 tor tne ha 



o mass, where r34o is 



the radius of the halo where the local density is p = 340p. 

The condition for screening of an isolated halo (i.e. 
not taking the environment into account) can be found 
theoretically. For a spherical top-hat over density with 
radius r34o and density p = 340p,„Q we have 



Am(?"34o) 
A Max 



1 



a > 1 

a < 1 



(33) 



where a is given in Eq. ( 10 ). The condition for screening, 
a > 1, can be written as 



M: 



340 



where ^340 = 



is the halo mass. This con- 



^340 = ^^^-340 

dition is not accurate for rea l halos as the top-ha t 



approximation is very crude (Clampitt et al. 2011) 
but it nevertheless is able to capture the essence ot the 
screening mechanism. In Table (121) we show the ratio 



of halos more massive than 5 x lO^^Mgun/ft, which are 
expected to be screened in the different simulations 
based on this simple approximation. The effectiveness 
of the symmetron screening mechanism increases with 
increasing halo mass il/340 and decreasing symmetry- 
breaking redshift zssb- This is because a larger value of 
zsSB means the symmetry is broken at higher matter 
densities and consequently a larger halo mass is required 
to restore the symmetry. 

In Fig. ([2]) we show Am (''340) for our simulations 
as a function of the environment for both large (blue 
circles) and small (purple circles) halos. Firstly, we note 
that the predictions from simulations with different /3 
are very similar and the only real effect of changing /3 
is in changing the normalization factor A^j^^ = 2fP . 
This can be understood from noting tha t ch anging /3 
only affects the geodesic equation Eq. (15) and not 
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the Klein-Gordon equatioij^ Eq. ([3|. However, for 
simulations with larger /3 we will on average have more 
massive halos because the fifth-force, and therefore the 
matter clustering, is stronger. This effect, which for our 
simulations is very small, can also be seen in Tabic (pi). 
Secondly, we note that the predictions of Am(?'34o) for 
high-mass and low-mass halos in low-density environ- 
ments are easy to separate at la for all our simulations. 
The small dispersion about the solid curves (which 
represent the average values in the two mass bands) seen 
in Fig. (I2]) is due to the difference in the halo masses 
within each mass band. 

To see this more closely, we have used the size of the 
circles (which represent halos) to denote their masses: 
bigger circles are more massive halos. We can see the 
clear trend that Am decreases with increasing circle size 
(or halo mass), and this confirms that in a given envi- 
ronment the screening of a halo, or equivalently Am, de- 
pends very sensitively on the mass of that halo. Thirdly, 
for very high density environments D — >■ we recover GR 
independent of the halo mass for all of our simulations, 
which is because the local value of (j) in a very high den- 
sity environment (which is often a part of or very close 
to a very massive halo) is small and so the fifth force is 
suppressed. 

In Fig. (|3) we show Am(?'34o) for our simulations as 
a function of the halo mass in both high-density (purple 
circles) and low-density (blue circles) environments. This 
figure shows that GR is recovered for larger halos, inde- 
pendent of the environment, as expected from Eq. (34). 
For low-mass halos we see a significant dispersion of ZIa/ 
from to the maximum value obtained in low-density 
environments for the same mass-ranges. This is because 
low-mass halos cannot efficiently screen themselves and 
must rely on the environment to get the screening. The 
environment, defined by D, ranges from i3 = up to 
D = 10 for these halos, and the lower the value of D the 
better screened the halo will be. To observe this point 
more clearly, in the figure we have also used the size of 
the circles to denote the value of D: the bigger circles are 
halos in environments with larger D (or lower density) 
and vice versa. As expected, we see a clear trend that the 
small halos which are efficiently screened geneally reside 
in high-density environments, while those which are less 
screened lie in low-density environments. 

Massive halos on the other hand can screen themselves 
efficiently and the environment only plays a small role in 
their total screening. 

In Fig. ^ we show AAf(r) as a function of the dis- 
tance r from the halo centre, for small and large halos in 
high and low density environments respectively. Again 
we see a large difference between large halos in dense en- 
vironments and small halos in low density environments. 
The r-dependence of Am {r) is seen to be rather weak in 
high density environments since the value of the scalar 
field inside the halo is mainly determined by the envi- 
ronment, while in low-density environments the value 



1 

2aI 



The Klein-Gordon equation Eq. (jsjl can be written as Di/i 



LPmo{l + ZSSB) 



^ — 1 -I- i/'^j where -0 = <l>/4'i 



iQ. Using this vari- 
2 

?/)VV'. Thus the 



able the fifth-force can be written F = f3 I - 

\ - -pi / 

solution ip is independent of /3 and its only the geodesic equation, 
through the fifth-force, which has a /3 dependence. 



of the scalar field mainly depends on the mass of the 
halo, which leads to a stronger r-dependence. Note also 
that in all the figures above the deviation from GR is 
stronger for higher sym metry-breaking redshift ^ssb , as 
expected from Eq. ( 34 ) , and for larger values of the cou- 



pling (3, which implies a stronger fifth-force and therefore 
a stronger effect. 

It should be emphasized that the environment depen- 
dence which are seen in the figures above will depend 
on the way the halos are binned, i.e. our definitions in 
Sec. (|4f. 

The results shown above are for halos at redshift z = 0. 
Another signature in the symmetron scenario is the red- 
shift dependence of Am- For halos at large redshifts, 
z > zssB, we have Aj\/ w 0, independent of the envi- 
ronment and mass as the symmetry has not been broken 
at the background level and cj) ^ almost everywhere 
in spac e, as demonstrated in Fig. (7-10) in Davis et al 



( |2011a[ ) . If this signature is found in observations, then 
dividing the observational samples in bins according to 
redshift one can probe the value of the symmetry break- 
ing redshift ^ssb • The maximal strength of the deviation 
will again probe /3, which can help to distinguish the 
symmetron from f{R) gravity. If one for example finds 
A^f-^ > I then f{R) cannot account for the deviation. 

6. COMPARISON WITH ANALYTICAL RESULTS 
Dur i ng the completion of this work a paper [Clampitt] 



et al. (20111) came out with semi-analytical halo scale 
predictions ifor the symmetron. In their analysis they 
assumed a NFW profile and calculated the symmetron 
fifth-force for isolated halos. The quantity of inte rest 



is (see Clampitt et al. (2011); Schmidt (2010) for 



details) which in our notation is given by 



ffv 



1 



/ r^p{r)FNAM{r)dr 
J p{r)F]<^dr 



(35) 



where Fjv is the gravitational force, g^-^j. is the average 
force to the average gravitational force over the halo. 
Since galaxies are spread around inside the halo, a mea- 
surement of the velocity dispersion of galaxies would 
therefore measure such an ave rage. One of the cases 
shown in Clampitt et al. (2011 ) can be compared to our 
simulation results, and as a consistency check we perform 
this comparison. 
In Fig, ([s]) we show together with the predictions 



from Clampitt et al. (2011 ). The results from their anal- 
ysis seem to be in good agreement with our numerical 
results. It would be interesting to see if their analysis 
can be extended by taking the environment into account 
to obtain the simulation results we have presented here. 
This would allow for an easier comparison with future 
observations as N-body simulations are in general very 
time consuming. 

7. SUMMARY AND CONCLUSIONS 

We have studied the environment dependence of the 
masses of dark matter halos in the symmetron modified 
gravity scenario. The potential governing the dynamics 
of the matter fields ($_ + $-|_) can differ significantly 
from the lensing potential <f>+ in this model, which leads 
to a clear difference between the mass of the halo as ob- 
tained from dynamical measurements and that obtained 
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Parameter: 


Value: 




0.267 


Simulation: 


A 


B 


c 


D 


E 


F 




0.733 


^SSB 


0.5 


0.5 


1.0 


1.0 


2.0 


2.0 


h 


0.719 


/9 


0.5 


1.0 


0.5 


1.0 


0.5 


1.0 




0.801 


L 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


Us 


0.963 




-^particles 


256^ 


Bq 


64Mpc/h 



TABLE 1 

The symmetron parameters for each of the N-body simulations A-F together with some re levant cosmologica l- and 

N-BODY parameters (SAME FOR ALL SIMULATIONS). FOR A COMPLETE LIST OF PARAMETERS SEe |DAVIS ET AL.| ( |2011a[ | . 



Simulation: 


A 


B 


C 


D 


E 


F 


Screened: 


99.82% 


99.83% 


64.38% 


65.27 % 


13.91 % 


14.12 % 



TABLE 2 

The percentage of halos more massive than 5 x IQ^^ M^ „^ I h which are expected to be screened (to some degree) in the 

DIFFERENT SIMULATIONS BASED ON THE APPROXIMATION EQ. (|34[| . THERE IS A VERY SMALL DIFFERENCE BETWEEN THE SIMULATIONS 
WHERE P = 0.5 (A,C,E) COMPARED TO /3 = 1.0 (B,D,F) EVEN THOUGH THE SCREENING IN ONLY SENSITIVE TO L AND ZsSB- THIS SMALL 
DIFFERENCE COMES FROM THE FACT THAT SIMULATIONS WITH STRONGER (3 WILL IN GENERAL HAVE MORE MASSIVE HALOS. 




• Halo mass M|_, log^^Q D> ^ • Clustered halos, • Isolated halos, 

# Halo mass M|_, log^^Q D< 1 Io9l0 ^ I°9l0 ^ 

Fig. 1. — The 3D distribution of the halos in the simulation box for a simulation with ^SSB = 0.5 and fi = L = 1.0. Left: The blue 

(D < 10) and purple {D > 10) spheres are 500 randomly selected halos in the mass range 11.5 < log^g j^.^^^ jy^ < 12.5. Because of the 

definition of the environment Eq. l |32[ l the largest halos in the simulation will almost always be in a low-density environment and we 
therefore separate our the 50 most massive halos and show these separately in black with a fixed size for the spheres. Middle and right: 
The value of Ajv/ for the clustered halos (logjQ D < 1) and the isolated halos (log^g D > 1) respectively. 



from gravitational lensing. Such an effect found in the 
symmetron model can be significantly stronger than in 
/(i?) gravity. This signature, which is unique to modified 
gravity, can in practice be measured by combining dy- 
namical (e.g. velocity dispersion) and lensing mass mea- 
surements of clusters of galaxies or even single galaxies. 
We find that the environmental dependence is strongest 
for small halos as very large halos are sufficiently massive 
to be able to s creen themselves. This im plies that using 
dwarf galaxies ( Jain fc Vanderplas|2011 1 might prove the 
best way to probe this ettect. 

This discovered feature of environmental dependence 
also allows us, in principle, to distinguish between dif- 
ferent modified gravity scenarios such as f{R), more 
general chameleons, DGP and the symmetron. In both 
DGP and f{R) the maximum fraction of the fifth- force 



to the Newtonian force in halos are around 30% while 
in chameleon/symmetron scenarios this fraction can be 
either smaller or larger, depending on the value of the 
couphng strength (3. DGP differs from f{R) and the 
symmetron in that there is basically no environmental 
dependence. There is also the possibility of measuring 
the redshift evolution of this effect by measuring clusters 
at high and low redshifts. As the symmetron force is neg- 
ligible for z > zggB we will recover the GR predictions 
for all clusters, independent of mass and environment, at 
high redshifts. 

Since different modified gravity theories can be highly 
degenerate with regards to both background cosmology 
and the growth rate of linear perturbations, it is crucial 
to identify new probes which can be used to separate 
them from each other. If one of these models are real- 
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zSSB = 0.5 



X 



P= 1.0 




(3= 0.5 



Log 10 D 



zSSB = 1.0 



Log-Lo D 



P= 1.0 




|5= 0.5 



Log^Q D 



zSSB = 2.0 



Log^o D 




Log^o D 



Log 10 D 



Fig. 2. — Am(''34o)/A5^'^'' for as function of the environment for large halos { jj^'"/^ > 5- 10^^, blue) and small halos {jj^^^ < 2-10^^, 

purple) where A^"'' = 2/3^ and where the size of the circles increases with the mass of the halos. The error bars are la. We see a clear 
difference between the value of Ajv/ between what we have defined as high-mass and low-mass halos. 
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zSSB = 0.5 



X 



o 



P= 1.0 



p= 0.5 




sun 



zSSB = 1.0 



sun 



o 



Log]_o M/M 



zSSB = 2.0 



Log^O M/M 

sun 





ID 



o 



P= 0.5 



Log-|_Q M/M 



Log]_QM/M 



Fig. 3. — AM{f:i4,o)/ ^^/^^ as a function of the halo lensing mass for high density environments (D < 10, purple) and low-density 
environments {D > 10, blue) where A^''^ = 2/3-^ and where the size of the circles increases with D (i.e. a smaller circle indicates a denser 
environment). The error bars are Icr. For the high-mass halos we recover GR independent of the environment as the effectiveness of the 
screening increases with mass. 
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Fig. 4. — AM{r)/A^j'^^ as a function of the rescaled halo radius r/r^io, where A^''" = 2/3^, for high (D < 10) and low (D > 10) 
density environments and small (red) and large (blue) halos. The error bars are Icr. For comparison we show the profile for all halos in the 
simulation (dashed green) in each plot and this curve has been displaced —5% in the r-direction to more clearly see the error bars. 
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12.0 12.5 13,0 1i6 n.<t 



Log^^O M/M 

Fig. 5. — g^;^ as function of halo mass for the isolated halos 
{D > 10) in a simulation with {/3 = 1.0, L = 1.0. zrrr = 0.5| com - 
pared with the semi-analytical results of | |Clampitt et al.] | |2011[ l, 
Fig. 3) (red). The dashed black line shows the UK prediction 
g^^^ = 1. The semi-analytical results are for isolated halos which 
agrees very well to the maximum g^;^ i n our simulations . Not e 
that the symmetron parameters used in [Clampitt et al.] |20TTJ: 
{/3 = 1.0, L = 1.2, zsSB = 0.54} are slightly ditterent compared 
to our simulation, and also for the definition for the virial mass 
M = M300 was used apposed to our M = M340. 



ized in nature then only a combination of many different 
probes will be able to pin down the correct theory. How- 
ever, the first step would be to detect a deviation from 
GR and a detection of the effect considered in this paper 
will be a smoking gun for modified gravity. It will there- 
fore be very interesting to look for this effect using data 
from upcoming large-scale structure surveys. 
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